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DETECTION  OF  SIGNALS  IN  NON-GAUSSIAN  CORRELATED 
NOISE  DERIVED  FROM  CAUCHY  PROCESSES 


INTRODUCTION 

The  detection  of  additive  signals  in  correlated  noise  or  radar  clutter  is  a  problem  that  is  of 
interest  and  has  been  extensively  studied.  The  optimal  detector  when  the  correlated  noise  is 
Gaussian-distributed  is  the  Wiener  filter,  or  matched  filter,  after  prewhitening  followed  by  a  thresh¬ 
old.  This  detector  can  be  obtained  by  using  the  Neyman-Pearson  procedure  that  maximizes  the  prob¬ 
ability  of  detection  for  a  given  probability  of  false  alarm  for  a  binary  hypothesis.  In  applying  this 
procedure  to  non-Gaussian,  correlated  noise,  three  problems  are  encountered.  First,  we  seldom  know 
or  can  easily  measure  the  required  multivariate  probability  density  of  the  noise;  second,  often  there 
are  unknown  parameters  that  must  be  accounted  for  in  some  way;  and  third,  the  likelihood  ratio 
obtained  in  the  test  sometimes  is  difficult  to  simplify.  All  three  of  these  problems  are  addressed  in 
this  study  for  a  given  situation  to  be  described. 

The  most  difficult  problem  encountered  is  obtaining  the  multivariate  probability  density  of  the 
noise.  A  procedure  for  constructing  an  approximate  representation  of  the  multivariate  probability 
density  is  described  by  Martinez,  Swaszek,  and  Thomas  [1].  The  procedure  constructs  the  desired 
multivariate  density  from  one  that  can  be  analytically  represented,  such  as  a  Gaussian  one,  by  using  a 
nonlinear  transform  to  map  the  one  into  the  other.  The  mapping  is  adjusted  so  that  the  marginal  dis¬ 
tributions  and  the  first  two  moments  of  the  constructed  multivariate  distributions  are  correct.  Often 
these  are  the  only  properties  of  the  clutter  that  can  be  measured  easily.  In  Ref.  2,  these  results  are 
extended  to  include  complex  numbers  to  represent  radar  baseband  signals,  to  provide  a  suitable 
transformation,  and  to  give  the  closed-form,  multivariate  probability  density  for  both  correlated 
Weibull  and  log-normal  clutter  that  is  correct  in  the  marginals  and  covariance  matrices.  The  purpose 
of  this  report  is  to  repeat  the  procedure  outlined  for  Ref.  2  except  using  an  underlying  Cauchy  pro¬ 
cess  rather  than  an  underlying  Gaussian  process.  The  examples  in  Ref.  2  use  a  bivariate  distribution, 
and  consequently  that  is  what  is  developed  here.  The  two  processes  are  identical  in  the  marginals  and 
the  first  two  moments  but  differ  in  other  respects. 

Once  the  bivariate  distributions  are  known,  the  optimal  detector  for  additive  signal  can  be  found 
in  closed  form  by  using  the  Neyman-Pearson  test.  In  Ref.  2,  an  approximation  to  the  Neyman- 
Pearson  test  was  found  that  simplified  the  test  and  removed  the  need  for  knowing  the  signal  amplitude 
and  phase.  This  detector  was  shown  to  perform  better  than  the  classical  matched  filter  after  pre¬ 
whitening,  which  is  optimal  only  for  Gaussian  noise. 

We  describe  a  means  of  generating  correlated  Cauchy  distributed  samples,  then  develop  an 
expression  for  the  bivariate  density  function  for  the  Cauchy  process.  This  is  followed  by  the  develop¬ 
ment  of  the  Weibull  and  log-normal  bivariate  probability  densities.  The  Neyman-Pearson  test  is  con¬ 
structed  and  simplified,  and  an  example  is  given.  Further  work  will  involve  comparing  detection 
results  of  the  detectors  described  in  Ref.  2,  and  those  developed  here  by  using  each  other’s  noise 
models. 


Manuscript  approved  August  8,  1987. 
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CONSTRUCTING  THE  BIVARIATE  CAUCHY  PROCESS 

Given  that  one  can  easily  construct  independent  Gaussian  random  variables,  Cauchy  variables 
can  be  constructed  from  them  through  a  nonlinear  transformation 

Zrc(l)  =  h~l(Zg(l)) 

and 


zlc(l)  =  h  l(z‘g(l)), 

where  zrc{l),  Zlc(l),  zg(l),  and  zg{l)  are  the  inphase  and  the  quadrature  components  for  the  Cauchy  and 
Gaussian  random  variables  respectively,  is  the  nonlinear  function,  and  /  indicates  the  / th  sam¬ 

ple.  The  inverse  transformation  h{-)  is  defined  by 

zrg(D  =  h{zrc{l) ) 

and 


z‘(l)  =  h(z‘c(l)). 

The  nonlinear  functions  are  specified  by  the  requirement  that  the  distribution  functions  over  each  vari¬ 
able  must  be  equal  and  are 

MzA /))  z  (/) 

pzg  (zg(l))dzg(l)  =  (zc(l))dzc(l) 

and 


z  (/»  h-'(zg(D) 

J '  PMDXk'U)  =  f 


The  Gaussian  density  is 


PMD)  =  7^—  exp 


Tff. 


2  a} 


and  the  Cauchy  density  is 


Pz(zc(D)  = 


where  a  is  the  standard  deviation  of  the  Gaussian  noise  and  b  is  the  Cauchy  parameter.  After 

Z8 

integration  the  nonlinear  function  is 


h  1  (*/(/))  =  b  tan 


} 


h  {zc\l )) 


'p2  erfc  1 


} 


2 


NRL  REPORT  9086 


where  the  complementary  error  functions  are  defined  as 


2  °°  2 

erfc(x)  =  —r — ■  (  e~'  dt 

Vtt  jx 

and  the  dot  represents  the  in  phase  r  and  quadrature  i  components.  The  integrations  are  classic  ones 
in  that  the  integration  over  the  Gaussian  density  yields  an  error  function  and  over  the  Cauchy  yields 
an  inverse  tangent  function. 


The  next  step  is  to  correlate  the  data.  Since  the  sum  of  Cauchy  processes  is  still  Cauchy  (shown 
subsequently),  the  data  are  correlated  by  using  some  of  the  same  samples  in  sums  for  each  of  the  new 
random  variables.  Furthermore,  the  results  are  to  be  obtained  for  both  the  independent  inphase  and 
quadrature  channels.  The  following  sum  operations  over  K  variables  accomplish  these  goals: 


rr 


1  K 
Kl= 1 


1  £  + 1 

Xr2  =  Y  L  Zrc(l), 
K  1=2 


xl 


and 

K  1  =  1 


1  +  1 

4(0- 

K  1  =  2 

The  marginal  probability  densities  of  each  of  these  new  random  variables  x\ ,  xr2 ,  x\ ,  and  x‘2  are  all 
identical  and  are  Cauchy,  given  by  the  density  function 


The  bivariate  density 


;  (xri  ,  X2  ,  x\  ,  x‘2)  —  Pxrxr  (x\ ,  X2)  >  X1  (x\  ,  X2) 


is  derived  in  the  next  section.  The  first  samples  x\  and  x\  can  be  made  highly  correlated  to  the 
second  set  of  samples  x2  and  x‘2  if  K  is  large.  Note  that  the  moments  are  not  defined  for  the  Cauchy 
process. 


BIVARIATE  CAUCHY  PROBABILITY  DENSITY 

The  bivariate  Cauchy  probability  density  is  found  by  using  characteristic  functions.  Because 
Pr  xr  (xj ,  x2)  and  Px, i  (x) ,  X2)  are  identical  in  form,  only  the  one  for  the  inphase  component  is 

developed.  The  characteristic  function  of  {x\,x2)  is  first  found,  then  Px^xr  (x^x^)  is 

found  by  using  the  inverse  Fourier  transform  of  it. 
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The  characteristic  function  Cxr  xr  (jw , ,  jw2)  of  Pxr  xr  (xj ,  x2)  is  defined  by 


where  E  is  the  expected  value  with  respect  to  x\  and  xr2 .  Substituting  for  x\  and  xr2  ,  the  characteris¬ 
tic  function  becomes 


Or,  jw2)  =  E  I  e 


Wj  K  w2  ^  +  1 

“F  E  Zed)  j~r  E  zrc{i) 

A  l  ~  1  e  A  l  =  2 


where  now  the  expected  value  is  over  zrc{l )  for  /  =  1, - K  +  1.  This  expression  can  then  be  writ¬ 

ten  as 


Cx\,x\  ObM)  =  C7rA 


z'cQ) 


jw_  l 

/if 


S&a:  + 1) 


jw 2 

a: 


/?2 


;(wi  +  w2) 

K 


by  noting  that  all  the  z£(/)’ s  are  independent  and  that  each  of  the  (K  +  1)  terms  are  by  definition  a 
characteristic  function.  The  characteristic  functions  (•)  are  characteristic  functions  of  Cauchy 

processes  and  are 


Czrc{l) 


w  | 


w, 


=  exp 


K 


CZrc(K  +  1) 


=  exp 


w2  |  b 
K 


) 


CZrc(D 


j(wx  +  w2) 


=  exp  i  — 


wl+w2\b 

K 


The  desired  characteristic  function  then  is 


Cxrxr  CM,  yw2)  =  exp 


+  (K  —  1)  |  W\  +  w2 1 ) 


The  bivariate  probability  density  is  obtained  by  taking  the  inverse  Fourier  transform  of  the 
characteristic  function  and  is  defined  by 


Px\,xr2  =  -^2  i  1  Cx\,x'2  (jw\,  jw2)  exp  j  -j{wlX\  +w2xr2)  idwldw2. 
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The  integration  is  performed  in  parts  over  different  regions  because  of  the  absolute  value  signs  in  the 
characteristic  function.  The  regions  are  defined  by  the  boundaries 


70  :  Wj  >  0, 

w2  >  0 

7'i  :  Wi  <  0, 

w2  >  0, 

|wi|  >  \  w2\ 

7' 0  :  Wi  <  0, 

Wj  <  0 

72  :  Wi  >  0, 

w2  <  0, 

|wi|  <  |w2| 

7j  :  Wj  >  0, 

w2<0,  |wj|  >  \  w2\ 

7'2  :  W!  <  0, 

n>2  >  0, 

|wi|  <  |w2| 

and  70,  7'0,  I\,  /' i,  7 2  and  7'2  are  the  values  of  the  integrals  over  those  regions.  Consequently 

Px'uxr2  (xl’x2)  =Io  +  I'o  +  I\+I'l+I2  +  I'2- 

The  integrations  are  sketched  out  in  Appendix  A.  The  result  is 


b 


fl 

,b  , 


1  + 


x\ 


1  + 


x 2_ 


(l  ~  Px )  ~ 


r 

ci 


x 2 


(1  ~  Px)2  + 


r  1 2 

,b  , 


1  + 


x\ 


(1  -Px)~ 


V b  , 


/""  "N 


ZL 

b 


(1  ~  Px)  + 


v_ b  J 


1  + 


x 2 


V.  J  J 


+ 


Pjt  ( 1  Px  ) 


"A 

*2 


x\  xr2 


f  ~A 

2  " 

2  " 

x2 

x\ 

*2 

(l  -  Px)2  + 

.b  J 

P  2  + 

l* 

Z> 

+ 


Px ( 1  Px ) 


rx2  ^ 


JC2 


(i  _  p*)2  + 

where  px  =  2/K.  A  small  <5x  indicates  a  large  correlation. 


r  -a 

2  " 

f 

-N 

2 

1— 

2  , 

x\ 

x2 

U  j 

Px  + 

l  b 

'  b  , 

(1) 
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The  density  ^  (x\ ,  jj)  is  identical  to  (xj ,  xr2)  except  superscript  r  is  replaced  with  i. 

The  density  was  checked  to  see  if  a  simple  Cauchy  was  obtained  after  numerical  integration  over  one 
of  the  variables.  This  was  obtained  for  a  number  of  cases.  This  fact  is  also  sufficient  to  guarantee 
that  the  area  under  the  density  is  one.  These  numerical  integrations  served  as  checks  to  the 
mathematics. 


BIVARIATE  WEIBULL  AND  LOG-NORMAL  PROBABILITY  DENSITIES 

Two  noise  distributions  of  interest  are  specified  by  the  Weibull  or  log-normal  marginal  distribu¬ 
tions  on  the  magnitude  of  the  signal.  To  map  the  bivariate  Cauchy  variables  just  derived  to  these  dis¬ 
tributions  on  the  marginals,  the  distribution  of  the  magnitudes  of  the  random  variables  given  by 

\xk  |  =  V  (x[)2  +  (xlk)2 

must  first  be  found.  Then  the  nonlinear  mapping  must  be  specified. 

The  distribution  of  the  magnitude  and  phase  of  the  Cauchy  random  variables  is  given  by 

P  i*  | ek  ( I  xk  |  ,6k)  =  JXt  Pxr(x[) Pxi  (4) , 
where  Jr  is  the  Jacobian  transformations.  Given 

4  =  \xk  I  cos  °k 


4  =  I  xk  I  sin  ek  , 


and 


dxk 

H 

d\xk  | 

d\xk  | 

Jx  = 

xk 

dxk 

dxlk 

ddk 

ddk 

= 


f*  % 

2 

7 xb 

1  + 

x{ 

u 

_ 

L  b  J 

**«)  = 


irb 


1  + 


4_ 


then 


^\xk\,ek  (lx* 


h)  = 


1  Xk  I 
7 r2b2 


1 

4 

1  + 

I-**  1 

2 

+ 

r 

H 

4 _ 

4 

sin2  2dk 

b 

V2  b  , 

V_  J 
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Integrating  over  6,  from  —  ir  to  tt  yields 


(I**  I ) 


where  the  integral 


4I  xk  1 

r 

1 

7T 

■5 

{\xk\2+l  b2)\j  1  + 

r 

\Xk  1 
b 

J 

> 

2 

J 

f 


dx 


a  +  b  sin2x  ^ a  (a  +  b ) 


tan 


-l 


r  _  -n 

^  /  a+b 

- tanx 

n 


found  on  page  152  in  Ref.  3  is  used. 


The  mapping  used  to  map  the  bivariate  Cauchy  into  a  bivariate  Weibull  is 


r  8  (1**1 )  r 
y[  =  - ; - : —  xl 

\*k\ 


and  the  inverses  are 


i  8  (|**|)  , 

yl  =  - : - ; - xl 

U  I  * 


r  _  gXtyfcl)  r 


= 


|y*  I 


yk 


(2) 


4  = 


«(|y*|)  ; 


\yk 


yic 


(3) 


for  k  —  1  and  k  =  2.  The  nonlinear  functions  are  specified  from  the  requirement  that  the  distribu¬ 
tion  function  on  |y^|  and  \xk  |  must  be  equal.  This  requires 

.s<|y*l)  ,  ,\yt\ 

J0  ^|*|(l**K<*h*l  “J0  P\yk\(\yk\)d\yk\ 

and 


J0  =  J0  p\yt\(\yk\)d\yk\ . 


The  Weibull  distribution  is  given  by 

=«  in  2 


I  yk 


a  —  1 


exp 


In  2 


I  yk 


M,, 
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where  My  is  the  median  value  and  a:  is  the  Weibull  parameter  that  ranges  from  —  .5  to  2  for  radar 
clutter.  For  a  —  2,  the  Weibull  reduces  to  the  Rayleigh  distribution.  The  lognormal  is  given  by 


P\yk  |(|y*l)  = 


<*i  m 


exP  — 


1 

2  In 

r 

Uk  1 

2 

2  of 

My 

V-  y  J 

> 

J 

where  My  is  the  median  value  and  07  is  the  standard  deviation  of  (ln|y*  |)2.  By  performing  the 
integrals  over  the  defined  marginal  densities,  the  nonlinear  functions  are  found  to  be 


8(\yk  I) 


8~\ 1**1)  =My 

for  the  Weibull  distribution  and 


In  •< 


2 

1  -  —  tan-1 

V  1 JC*  I2  +  b2 

7T 

l  b  J 

- 

ln2 


1/a 


g  ( I  yk  I )  =  b 


■4 


tan  —  erfc 
8 


V2 

°l 


In 


I-** 


M, 


-  1 


8  (|**  |)  =  My  exp  " 


gf 

V2 


erfc  1 4 


1  -  Ltan_1 

7T 


'Tm  2  +  ^2 


for  the  lognormal.  The  function  erfc  is  the  complimentary  error  function,  erfc  1  is  its  inverse,  and  In 
is  the  natural  logarithm.  Appendix  B  provides  the  details  for  finding  these  functions. 

Finally  the  new  bivariate  Weibull  or  bivariate  lognormal  density  can  be  specified  in  terms  of  the 
bivariate  Cauchy  by 


Py\,y'2,y\ ,  y‘2  0>1  ,  /l  ,  *2 .  ?2>  =  Px'.xi  (*1  »  xt)PjA  .  xk  C*‘l .  *2  VyCfl  ,  ?  1  ,  ?2  »  )£) 


x\ ,  x'2 


x\ ,  x2 


g(bil) 

*\  =  \>,\  y\ 


x{  =  g(|yt|)  y[ 

Ini 


(4) 
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where 


Jy(yri,  y\,  yr2,  yli)  = 


o 

j2 


=  j,j 


\J2 


and 


dxl  dxl 


dy[  dy[ 


dx[  dxk 


g(|y*  l)g'(|yt  I) 

|y*  I 


dy'k  dyl 

where  g  '(•)  is  the  derivature  of  g(-)  with  respect  to  \yk  \ .  The  Jacobian  for  the  Weibull  is 


Jh  = 


ob2  ln2 

n»n 

a  —2 

r  “> 

7T 

r  -\ 

7 r 

r 

^ _ 

VMy  J 

i/'tan 

- d 

2 

J 

sec 

—  ~yb 

2 

^  J 

where 


7 T 

T 


—  In2 

1  y*  1  /My 

i 

L 

L  J 

j 

and  for  lognormal 


-&2V 5t 

I2 


1  ^  , 

exp  —  •<  - In 


|y* 


where  \p  for  this  case  is 


7T 


t  =  ~  erfc 


V2 

<*/ 


In 


AL 


|y* 


tan  sec2  \p 


AL 


(5) 


(6) 


The  bivariate  Weibull  distribution  was  checked  by  numerical  integration.  The  signals  were  con¬ 
verted  to  magnitude  and  phase  through  another  transform  and  then  integrated  over  the  two  phases  and 
one  magnitude.  The  remaining  marginal  density  was  Weibull  as  desired;  consequently  the  integration 
over  all  variables  is  one. 


NEYMAN-PEARSON  TEST 


The  classical  radar  solution  for  optimally  detecting  a  signal  in  additive  noise  is  the  Neyman- 
Pearson  procedure  for  a  binary  hypothesis  test.  This  processor  maximizes  the  probability  of  detection 
given  a  fixed  false  alarm  rate.  The  processor  is  specified  by  the  likelihood  ratio  and  in  this  case  is 


\ np 


Pyri,yr2,y\,y‘2  -sl)’  07 1  “  s[),  (yr2  -  s2),  (y2  -  sl2)) 


P.,r 


y\<y2>y\<yi 


(yi»  y  i  >  y{) 
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where 


s\  =  s  cos  (<t>s). 


■s'l  =  Ssin(<j!>s), 


sr2  =  S  cos  (cj>s  +  A)  ,and 


sl2  =  S  sin (<j)s  +  A). 

Here,  S  and  6S  are  the  signal’s  unknown  amplitude  and  phase,  and  A  is  the  phase  rotation  of  the  sig¬ 
nal  between  samples  1  and  2.  The  signal-to-noise  ratio  is  defined  to  be 

(S/AO  =  S2/ 2a]  , 

where 


a2  =  VAR  (y[)  =  VAR(y'k) 

and  VAR{')  is  the  variance.  A  very  lengthy  closed  form  expression  for  the  likelihood  ratio  can  be 
written  by  combining  Eqs.  (2)  through  (6).  However,  the  test  is  very  difficult  to  implement  and  con¬ 
tains  unknown  parameters.  Consequently  an  approximation  to  this  processor  is  next  considered. 

APPROXIMATE  NEYMAN-PEARSON  TEST 

The  noise  is  assumed  to  be  heavily  correlated  ( px  very  small),  and  the  signal  is  assumed  to  be 
small  with  respect  to  the  noise.  Equations  (2)  and  (3)  under  signal  condition  can  be  approximated 
from  Ref.  2  as 


r  sd?*!)  . 

xk  ~  — ; — ; — ~  csk 


:  *(|y*|)  ; 

XL  ~  — ; - ; — yi  —  csl . 

Ittl 


By  using  this  approximation,  the  likelihood  ratio  becomes 


px'r'Sxi » xi)  px\x^x\ -  4Vv(yi  y\ ,  ■  ■  ■ ) 


x\ 


s(bil)  r 
— i — r-y  \  ~csk 
yi 


^Np  ~ 


Px\,xrSX\  .  x2)Px‘,  .x‘-,(Xl  >  x2)Jy(y  1 .  yi.  ’  ’  ’) 


x\,x‘2 


x\ 


^(ItiI)  r 

b»l  Jl 
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Assuming  that  the  Jacobians  cancel  and  that  x\,x2,x\,  and  x\  are  the  transformed  samples,  the 
likelihood  function  is 


\ np 


Pxruxr2  [C*l  (X2  -CS$)]  Pxi,x>Vx\  -CS\),  (X^-Csi >)] 

Px\fx\  (Xl,  Xr2)Pxii  X,2(x[,  xl2) 


The  expressions  for  Pxr  xr  (x\  ,  x2)  and  xi  (x\ ,  x‘2 )  are  simplified  by  noting  that  the  factor 


Px  + 


2 


in  the  denominator  of  Eq.  (1)  that  defines  P xr  (x\ ,  x2)  basically  controls  the  central  peak  of  the 
probability  density,  since  px  is  small  and  x\  ~  x2  under  high  correlation  conditions.  Consequently, 


Px'ux'2  (*\’Xr2) 


7 


0  . 

r  -\ 

x\  x2 

2 

Px  + 

and 


Next,  \a  is  defined  by 


where 


Px\  jA  ’  *2> 


7 


7  1 

r 

H 

^ _ 

2 

Px  + 

lb  b  j 

Xa  XflXfl  ’ 


r  \2 


XI  = 


bz9lx+(x\  -xr2) 


b2Px  +  l(X\  -CS  i)  “  (xr2  -  CSr2 )] 


r  y.2 


and 


XI  =  _ bp*  +  (x i  -Xl2)2 _ 

b2pl  +  [(x:1!  -  cs\ )-(xi2-  cs^)]2 

These  can  be  approximated  by 


Xflr  =  1  +  2c (5 1  -sr2)(x[  - xr2 ), 
K  =  1  +  2c (s\  ~s\)(x\  ~x\). 

Then 


K  ~  (s'l  -4)C*1  ~*2)  +  (5l  -®2)W  -X2) 
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where  the  higher  order  term  has  been  ignored  and  the  other  coefficients  have  been  incorporated  into 
\a .  This  expression  states  that  the  data  are  transformed  back  to  the  bivariate  Cauchy  random  variable 
and  then  the  correlation  between  samples  is  removed  by  subtraction.  This  seems  reasonable  since  it 
was  in  the  Cauchy  domain  that  the  data  were  correlated.  Consequently,  the  correlation  probably 
should  be  removed  in  this  domain.  However  the  Cauchy  domain  is  not  an  easy  one  in  which  to  set 
the  thresholds  because  of  the  large  tails  on  the  distribution.  Since  the  Cauchy  random  variables  were 
generated  from  Gaussian  ones,  the  subtracted  random  variables  are  transformed  to  Gaussian  by  h . 
The  detector  finally  used  is  then  given  by 

K  =  -xr2)]2  +  mx\  -4)]2- 

The  detector  transforms  the  measurements  from  bivariate  Weibull  or  lognormal  to  bivariate  Cauchy, 
removes  the  correlation  by  subtraction  in  both  the  r  and  the  i  components,  transforms  the  remaining 
Cauchy  to  Gaussian,  and  computes  the  energy  in  both  r  and  i  signal  components.  This  is  intuitively 
satisfying  in  the  sense  that  the  processor  inverts  the  way  the  data  are  generated. 

SIMULATION  PERFORMANCE  OF  AN  EXAMPLE 

The  performance  of  the  detector  is  evaluated  by  using  data  samples  obtained  from  the  con¬ 
structed  multivariate  density  of  the  clutter.  The  performance  curves  are  defined  as  the  probability-of- 
detection  vs  signal-to-noise  ratio  given  a  fixed  probability  of  false  alarm  for  various  problem  condi¬ 
tions.  The  performance  of  the  approximate  Neyman-Pearson  test  is  compared  to  that  of  the  matched 
filter  defined  by 


Kf  =  |  S'R -‘Y I  , 

where  the  traditional  complex  notation  has  been  used  to  represent  the  data  vector 


y\  +  jy\ 
yr2  +  jy  2 


where  j  —  ^  —  l. 


The  covariance  matrix  is 


Ry  =  Oy 


1  Py 
Py  1 


where  py  is  the  correlation  coefficient  of  the  Weibull  noise  that  has  an  even  spectrum,  and  the  signal 
vector  is 


s\  +  js\ 

S  =  r  ,  .  i  . 

+  JS 2 

The  bar  represents  the  conjugate,  t  represents  the  transpose,  -1  represents  the  inverse,  and  |:| 
represents  the  absolute  value. 

For  this  example  N  =  2,  N  =  40,  b  =  1,  My  =  1,  aZg  =  1,  <t>s  is  arbitrary,  and  the  steering 
vector  is  sj  =  1/^2,  s\  =  0,  sr2  =  0,  and  s‘2  =  1/^2.  For  Weibull  clutter  of  parameter  7=1, 
1.5,  and  2,  it  was  found  by  simulation  that  ay  =  1.42,  .97,  and  .84,  respectively,  and  that  py  =  .972 
for  all  three  values. 
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The  data  used  in  the  simulation  are  generated  by  converting  independent  Gaussian  distributed 
samples  to  Cauchy  distributed  samples,  correlating  them  through  a  summer,  and  then  converting  them 
to  Weibull  distributed  samples  on  the  marginals.  This  data  generation  is  used  to  set  the  thresholds 
through  the  probability  of  false  alarm.  Signals  are  then  added  with  random  phase  to  the  bivariate 
Weibull  distributed  samples  to  obtain  the  data  required  for  generating  the  operating  characteristics  of 
the  detectors. 

The  probability  of  false  alarm  pja  vs  threshold  is  computed  by  using  Monte  Carlo  simulation 
techniques.  Random  samples  of  noise  y \  ,  y2 ,  y\ ,  yl2  are  generated.  A  normalized  histogram  of  the 
number  of  samples  having  values  in  a  small  interval  for  each  detector  is  found.  After  many  trials,  a 
curve  of  the  probability  density  of  the  output  of  each  detector  vs  Xn,  or  \mf  is  obtained.  The 
probability  of  false  alarm  is  defined  as  the  sum  of  all  values  of  the  density  from  X*,  or  \*mj  to  infin¬ 
ity.  For  small  values  of  the  probability  of  false  alarm,  an  importance  sampling  technique  is  used. 
This  procedure  distorts  the  generation  of  random  samples  so  that  more  false  alarms  occur  than  should 
and  then  compensates  for  this  in  the  weightings  used  to  generate  the  histogram.  This  technique  is 
outlined  in  Appendix  C.  The  threshold  settings  vs  probability  of  false  alarm  for  the  new  approximate 
detector  and  the  matched  filter  are  shown  in  Table  1  for  a  =  1,  1.5,  and  2.0  and  for  probability  of 
false  alarms  as  low  as  10  ~7.  The  reason  the  thresholds  for  the  approximate  detector  are  independent 
of  a  is  that  the  detector  maps  the  Weibull  distributed  measurements  to  an  identical  Cauchy  distribu¬ 
tion  regardless  of  what  a  is.  Although  the  thresholds  are  independent  of  a,  this  does  not  mean  that 
the  detector  performance  is  independent  of  a. 

The  performances  of  the  detectors  are  compared  by  observing  the  probability  of  detection  vs 
signal-to-noise  ratio  for  a  fixed  probability  of  false  alarm.  The  probability  of  detection  is  computed 
by  using  a  Monte  Carlo  simulation.  The  fraction  of  time  the  detector  output  exceeds  the  threshold  for 
a  set  of  randomly  generated  samples  is  computed.  This  number  is  the  probability  of  detection. 

Performace  results  are  shown  in  Figs.  1  to  3  for  Weibull  parameters  of  a  =  1.0,  1.5,  and  2.0, 
respectively.  In  all  cases,  the  probability  of  false  alarm  is  10-7.  Figures  1  to  3  show  that  the 
approximate  detector  performs  better  than  the  matched  filter.  In  all  cases,  the  new  approximate 
detector  performs  better  than  the  matched  filter;  that  is  optimum  for  stationary  Gaussian  bivariate  dis¬ 
tributions.  This  is  true  even  when  a  =  2  and  the  marginal  densities  become  Gaussianly  distributed. 
Even  though  the  marginals  are  Gaussian,  the  data  are  not  bivariate  Gaussian  and  consequently  the 
new  detector  is  better.  A  functional  flow  of  the  detector  is  shown  in  Fig.  4. 


Table  1  —  Probability  of  False  Alarm  vs  Thresholds  for  Weibull  Clutter 


Pfa 

a 

10-1 

10-2 

10~3 

10~4 

10-5 

10-6 

r- 

1 

O 

r-H 

1.0 

.51 

1.65 

2.52 

3.30 

3.85 

4.15 

4.35 

Approximate  filter 
X* 

1.5 

.51 

1.65 

2.52 

3.30 

3.85 

4.15 

4.35 

A a 

2.0 

.51 

1.65 

2.52 

3.30 

3.85 

4.15 

4.35 

Matched  filter 

1.0 

5.0 

30. 

75. 

115. 

155. 

195. 

230. 

x*. 

1.5 

9.0 

45. 

95. 

135. 

170. 

200. 

230. 

Amf 

2.0 

9.0 

54. 

104. 

137. 

165. 

189. 

210. 
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Fig.  1  —  Detector  operating  characteristics  for 
Weibull  clutter,  a  -  1.0  and  pfa  =  10-7 


Fig.  2  —  Detector  operating  characteristics  for 
Weibull  clutter,  a.  =  1.5  and  pfa  —  10~7 


Fig.  3  —  Detector  operating  characteristics  for 
Weibull  clutter,  a  =  2.0  and  pfa  =  10"7 
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Fig.  4 — Signal  processor  functional  diagram 


SUMMARY 

A  procedure  to  detect  a  target  in  non-Gaussian  correlated  noise  is  obtained.  Since  the  bivariate 
probability  density  is  unknown,  an  approximate  one  is  constructed.  The  constructed  density  matches 
the  true  density  in  the  marginals  and  first  two  moments.  The  mappings  required  are  found  for  both 
the  Weibull  and  the  lognormal  clutter  distributions.  The  bivariate  density  is  formed  from  first  map¬ 
ping  an  independent  Gaussian  random  variable  into  Cauchy  distributed  random  variables,  correlating 
with  summers,  and  then  mapping  so  that  the  new  bivariate  density  has  the  desired  marginals  and  first 
two  moments. 

A  Neyman-Pearson  test  was  obtained  for  detecting  an  additive  signal  in  this  bivariate  distributed 
noise.  Since  the  test  was  complicated  and  contained  signals  with  unknown  amplitude  and  phase,  it 
was  simplified.  The  simpler  test  transforms  the  correlated  Weibull  or  log-normal  measurements  back 
to  Cauchy  distributed,  removes  the  correlation,  and  finally  transforms  back  to  “Gaussian-like  distrib¬ 
uted”  variables  where  the  result  is  compared  to  a  threshold.  The  performance  of  this  detector  was 
compared  to  that  of  the  matched  filter  operating  on  the  same  data  set.  It  was  found  that  the  new 
detector  performed  better  than  the  matched  filter.  One  could  speculate  that  the  best  detector  always 
transforms  the  noisy  data  to  that  domain  where  the  correlation  is  best  removed  and  then  transforms  it 
so  that  stationary  Gaussian  random  variables  are  obtained  on  all  examples. 
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Appendix  A 

INTEGRATIONS  FOR  BIVARIATE  CAUCHY 


The  integrations  70  and  7'0  are 


h  = 


1 

(2tt)2 


oo  oo  — 

j  j  exp  {  —  b\wx  +  w2  +  (N  —  l)(wt  +  w2)] 


~jwixi  ~  jw\xi)  dwldw2 


1  o  o  — 

/'  0  =  - 7  (  exp[-b[-wl-w2-(N-l)(wl  +  w2)]-jwlxl-jw2x2}dwldw2 

(hr)2  J-°°  J-“ 

where  the  superscript  r  and  i  have  been  dropped  for  convenience.  By  changing  variables 
w\  =  —  vvj  and  w'2  =  —  vv2,  the  two  integrations  can  be  combined 


/o+/'o  = 


1 


(2tt)2 

This  can  be  simplified  to 


(  (  exp { —Nb (w i  + Wo) }  [exp{-j (x iwl+x2w2)}  +  exp[+j (x l+x2w2)}]dw ^dw2. 

Jo  Jo 


In  +  I'n  = 


2  f"  -Nbw,  ,  r  0 

L  e  COSJ^WjAv,  Jo 


(2tt) 


,2  .>0 


— Nbw 7  , 

e  cos  x2Widw2 


2  f 00  —Nbw{  .  }  f  °°  -Nbw 2  .  7 

- -I  e  smjc^^wij  e  sin  x2w2dw2 . 

(2ir)2  0  0 

These  integrals  are  standard  table  integrals  and  consequently 

_ (M>)2  —  x  xx2 _ 

[(Nbf+xl  ]  [(Nbf+xt] 


/n  +  /'n  = 


(2tt)2 


The  integrals  7 j  and  /']  are 


l  00  0  — 

7,  =  - 7  (  f  exp  {  -  b[wx  -  w2  +  (N  -  lXwj  +  w2)]  ~jwlxl  -jw2x2]  dw2dwx 

(2tt)2  j° 


1  c°  r._,Vl  “ 

I'l  =  - r  f  exp{-b[-wl  +  w2-(N-l)(wl  +  w2)]-jwlxx-jw2x2}dw2dw1. 

(2? r)2  J-°°  J° 

These  can  be  combined  after  changing  variables  to  yield 


|  OO  0  __  - 

/.+/'.=  - -f  f  exp{-Nbwx-(N -l)bw2}[txp{-j(xlwl+x2w2)} 

(2tt)2  J°  J_w> 


+  exp {+  j(xxw1  +x2w2) }  ]dw2dwx 
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This  can  be  simplified  to 


Ii+I'i  =  - 7 

(2tt)2 


oc 

Jo 


-Nbw, 


cos  wxyx 


1  +  (N  -  2)fcw2 


L  J0 


COS  W2y2dw2 


dwi 


2  » v“  —Nbw.  .  r 

H - r  \  e  sin  WiV i  1 

(2 it)2  J°  U° 

After  integrating  over  vv2,  I  x  + I'  x  become 


'  +(N  —  2)bw2  . 


sin  'w2y2dw2 


dwx . 


/i  +  /',  = 


(27r)2[(A-2)fc)2+x22] 


(A  —  2)fe  e  2fcvv‘  cos  —  x2)dwx 


-x^o 


—2  bw,  . 

e  sin  w 


i(xi  -x2)dwx  -  ( A  —2 )b  |o  e  Nbw'  cos  wxxxdwx 


+  x2  |Q  e  Nbw'  sin  wxxxdw 


After  integration  over  wx,  then 


/i+/'i  = 


_ 2  [2(A  —  2)A>2  —  -x2Qy !  --t2)] _ 

(2tt)2  { [(A  -  2 )b  ]2  +  xl  }  [(2b  f  +  (xx-  x2)2] 


2[N(N  -2)b2  -xxx2] 


(2tt)2  {[(A  -  2 )bf  +  xj }  ((Nb)2  +  ;c 2  ] 

The  final  integrations  for  /2  and  /'2  are  similar  and  are 

1  0  ~wi  — 

h  =  ■— jf  Loo  10  exp  i~b\-wi~w2-(N  ~  1)(”T  +w2)] -jWlXl  -  jw2x2 }  dwxdw2 
Y  oo  o  _ 

I’ 2  =  ;  2  In  f-w  ^P{-f’[-wl+w2  +  (N-l)(wx  +  w2)]-jwlxl-jw2x2]dwxdw2. 

(27r)  U  W2 

Combining  and  changing  variables  yield 


i2  +  r2  == 


(2tt)2 


OC 

Jo 


— Nbw , 


COS  WJJ2 


n 

J„ 


2  +  (N  —2)bW] 


cos  wxyxdwx 


dw0 


2  r°°  —Nbw  2  . 

H - t  \  e  sin 

Jo 


(2tt/ 


v 

W&2  [ 


2  +  {N  -2)bw,  .  , 

e  sin  wxyxdwx 


dwn 
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In  a  similar  manner 
tables;  the  result  is 


to  /j  +I'i,  first  one  integrates  over  wY  then  w2  by  using  standard  integration 


2[2(A^  -  2)fe2  -  jc  t(.x2  -  jc  t)] 
(2ir)2[[(N  -  2  )bf  +  x\}  [(2b)2  +  (x2-xx)2] 

2[N(N  -2)b2-xxx2] 

(2tt)2 {[(N  -2)b?+x\  [( Nb )2  +  x\  ]  ' 
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Appendix  B 

DISTRIBUTION  FUNCTIONS 


The  distribution  function  for  the  magnitude  of  the  Cauchy  variables  over  the  in-phase  and  quad¬ 
rature  components  is 


io 


4l*lt  1 

“\ 

7T  ^ 

(\xk\2  +  2b2)*J  1  + 

/T 

\xk  1 

& 

2 

► 

d  1**1  • 


Changing  variables  to 


£  =  \xk\2  +  b2 


yields 


,vi-,i,+‘i  l  4  , 

=  ,  dt  =  —  tan-1 

w  i  l2 


' +b  £*+b 


7 r 


vi 

1**  1 

|2  +  &2 

&  J 

-  1. 


The  distribution  function  for  Weibull  is 

a  —  1 


rtol  ,  J  l» 

f  a  In  2 
Jo 


and  is  obtained  by  first  noting 

M 

So  exp 


NL 


In  2 


exp 


r  ,  ,  -\ 


In  2 


I  yt 


Mv 


d \yk 


I  y* 


M, 


ln2 


|y* 


Then  we  obtain 

1  —  exp  -T  —  In  2 

The  distribution  function  for  the  lognormal  is 


y* 


My 

J  j 


I. 


l>*l 


I  y* 


exp—  -< 


1 


2ff;2 


2  In 


I** 


M, 


y  d  |  yk  |  , 
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which  is 


Jo 


^2^ 


iroi 


exp 


2a? 


2  In 


\yk 


Mv 


2  In 


I  yk 


v.  My 


This  is  the  form  of  an  integral  over  a  Gaussian,  which  is 


1  .  1  , 
—  +  —  erf 
2  2 


V2 


In 


I  yk 


Mv 


These  distribution  functions  are  used  to  define  the  mappings  g(-)  and  g  (■). 
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Appendix  C 

IMPORTANCE  SAMPLING 


The  importance-sampling  procedure  [C1,C2]  is  a  Monte  Carlo  technique  that  distorts  the  genera¬ 
tion  of  the  random  number  so  that  the  events  of  interest  occur  more  frequently  than,  but  in  the  same 
manner  as,  the  events  occur  in  nature.  The  probability  of  the  event  occurring  is  then  compensated  for 
with  a  weighting  factor  so  that  the  true  probability  of  the  event  is  obtained. 

To  compute  the  probability  density  of  the  detector,  the  filter  values  are  quantized  by 

X  =  m  AX  , 

where 


m  =0,...  M  -1. 

The  probability  density  is  computed  by 


1  N 

p(K  is  between  (m  —  1)  AX  and  m  AX)  =  —24 

™  k  -  l 

where  N  is  the  number  of  Monte  Carlo  samples  and  5*  =  1  if  no  importance  sampling  is  used.  The 
equation  is  simply  counting  the  percentage  of  time  the  samples  fall  in  the  mth  interval.  For  impor¬ 
tance  sampling  and  this  case,  the  first  and  last  sample  of  the  independent  Gaussian  variables  are  dis¬ 
torted  in  their  standard  deviation.  The  weight  factor  is  computed  by  ratioing  the  true  density  to  the 
distorted  density  evaluated  with  the  distorted  data  samples.  In  this  case 


°d 


exp 


_j_ 

2 


u;o)]2 + [(zid)]2 


bg(i)32  +  [4(i)]2 


exp 


[zUN  +  l)]2  +  [zl  (N  +  l)]2  [zrJN  +  l)]2  +  [z‘  (N  +  l)]2 


°d 


where  ad  _is  the  distortion_of  the  standard  deviation  on  the  first  and  last  samples.  The  values  z', ( 1 ) , 
Zg(l),  Zg(N  +  1),  and  z'(N  +  1)  are  the  sample  values  under  distortion.  This  distortion  decreases  the 
correlation  in  the  final  data  and  raises  its  variance  slightly.  Consequently  more  false  alarms  for  the 
same  threshold  are  produced  after  distortion  but  the  probability  is  reduced  because  of  8k . 
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